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The spherically symmetric magnetic monopole in an SU(2) gauge theory coupled to a massless 
Higgs field is shown to possess an infinite number of resonances or quasinormal modes. These modes 
are eigenfunctions of the isospin 1 perturbation equations with complex eigenvalues, En = uJn — iy-n, 
satisfying the outgoing radiation condition. For n — > oo, their frequencies uJn approach the mass of 
the vector boson, Mw, while their lifetimes l/"fn tend to infinity. The response of the monopole 
to an arbitrary initial perturbation is largely determined by these resonant modes, whose collective 
effect leads to the formation of a long living breather-like excitation characterized by pulsations with 
a frequency approaching Mw and with an amplitude decaying at late times as t~^^^. 



PACS numbers: 03.65.Ge, 11.27.+d, 14.80.Hv 

Magnetic monopoles - magnetically charged finite en- 
ergy solutions in field theories with spontaneously broken 
gauge symmetries [J] - play an important role in a num- 
ber of field theoretic considerations. They can account 
for the quantization of the electric charge, catalyze the 
decay of proton, might be among the relevant topological 
defects in the universe, determine supersymmetric vacua, 
etc. - see for reviews. 

In the present Letter we point out yet another in- 
teresting aspect of the monopole which seems to have 
gone unnoticed so far - the existence of its quasinormal 
modes (QNM) or resonant excitations. We show that 
the Bogomol'nyi-Prasad-Sommerfield (BPS) monopole 
admits an infinite number of QNM. They manifest them- 
selves as resonance peaks in the low energy scattering 
cross section of isospin 1 scalar particles. The QNM 
can also be described as complex energy solutions of the 
isospin 1 small fluctuation equations that are regular at 
the origin and satisfy the outgoing radiation condition at 
spatial infinity. For previous studies of the reai energy 
small fluctuations around the monopole we refer to Q . 

We shall demonstrate, in particular, that the QNM 
of the monopole lead to a universal late time behav- 
ior of the perturbed monopole by giving rise collectively 
to a quasiperiodic, long living excitation whose ampli- 
tude decays as t~^^^ at late times {t being the standard 
Minkowski time). In a recent numerical study of Fodor 
and Racz a breather-like excitation of the monopole re- 
taining a considerable fraction of the energy of the exter- 
nal perturbation has actually been observed It has 
been one of our aims to explain their observations. 

It is worth noting the striking analogy of the t~^^^ 
asymptotic behavior of the monopole with the late time 
evolution of massive fields in black hole spacetimes 0. 
Black holes also possess the QNM, which have been much 
studied recently (see e.g. '3|). 

We consider a Yang-Mills-Higgs (YMH) theory with 
gauge group SU(2) defined by the Lagrangian 
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field strength tensor, D^<^°- = 9^$" + SabcA^^"^ denotes 
the covariant derivative, and in our units the mass of the 
gauge bosons is equal to one. My/ = 1, while the mass of 
the Higgs particle is \/2A . The energy-momentum tensor 
of the theory defined by Eq.lQJ is 

We restrict our analysis to the 'minimal' spherically 
symmetric sector, where the ansatz for the YMH fields is 
given by = 0, 



A"=£a.fe^(l-mt,r)), ^- = '^H{t,r), (2) 

where a, i , A: = 1, 2, 3 and = x^x^ . With 0= 
the YMH equations reduce to 

{r'^n + W^ + - l)W = 0, 
{r'^D + 2W^ + X{H'^ -r'^))H = 0. (3) 

The static, finite energy solution of these equations is the 
't Hooft-Polyakov monopole 0. For the special case, A = 
0, the solution is analytically known - the BPS monopole, 
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H{r) = r coth ?' — 1 



(4) 



We shall consider small fluctuations around the static 
monopole background: W — > W{r) + w{t,r) and H 
H{r) + \/2h{t, r). Linearizing Eqs.© with respect to w 
and h, we obtain 
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- l)w = -q 2V2WH h, (5) 



\{3H^ ~ r^))h = -q2V2WH^ 



(6) 



jbc^u^!^ is the non-abelian 



Here an auxiliary parameter, q, has been introduced, 
presently q — 1. In this Letter we concentrate on the 
case A = 0, when W, H are given by Eq.l^J. 

Resonant scattering.— First we shall demonstrate 
that Eqs.llSJl,® describe resonance phenomena indeed. 
Separating the variables as w = 3?(e~'"*Wcj(?')) and 
h = 3fi(e~'"*/i^(r)) with real uo, Eqs.(|SJl,® become a 
standard two-channel Schrodinger system. As it will be 



FIG. 1: The scattering cross section a{ij). 



FIG. 2: The complex solutions of Eqs.Q,® for n = 3. 



clear from what follows, the frequency spectrum is con- 
tinuous, uj^ > 0. Regular solutions of Eqs.©,® have 
to satisfy the conditions ^ h^^ ^ for r ^ 0, and 
they can be normalized for r oo such that 

h^{r) ^ sm{Lor + Situ)) , w^{r) ^ CiLo)ri e-''\ (7) 

where v = Vl — uj'^- The h-&e\d is massless, so it oscil- 
lates as r cx) for any value of w. The w-field is massive, 
and for < 1 it shows a bound state type behavior, with 
exponential decay as r — s- oo. The fact that the it;- field 
is non-radiative for < 1 plays the crucial role in our 
analysis, and below we shall concentrate to this frequency 
range. Eqs.©,® describe in this case the scattering of 
a massless /i-radiation on the monopole surrounded by a 
confined massive ui-field. We note that for oj^ < 1 this is 
effectively a one-channel scattering problem, so that the 
scattering cross section is given by the standard formula 
ct(w) = {4TT/u;^)sm\Siij)). 

It is worth noting that the interaction of the /i-field 
with the monopole is in fact short range, so that the 
cross section is finite. We integrate Eqs.©,® numeri- 
cally to obtain Wc^{r) and hui{r) subject to the boundary 
conditions The resulting cross section a{uj) shown in 
Fig^exhibits a sequence of resonant peaks accumulating 
near the value lo = 1. This can be so interpreted that for 
certain energies of the incident /i-radiation the monopole 
core gets strongly excited. 

QNM — numerical results.— The scattering reso- 
nances can be usually related to the quasinormal modes 
- complex energy solutions satisfying the purely outgo- 
ing wave condition at r = oo. To construct the QNM, 
we integrate Eqs. Q , l|njl with w — 3fi(e~*^*W£;(r)) and 
h = 3fi(e~*^*/i£;(r)), where we and He are complex, the 
energy E = uj — i"f, and 

Ar^ ^ hE{r) ^ e*^^ Br^ ^ WEir) Crie-"'' (8) 

for ^ r — > oo. With the 'shooting to a fitting point' 
numerical method we find a discrete family of global 
solutions wsir), hsir) subject to the boundary condi- 
tions lO labeled by n = 1, 2, . . ., the number of nodes of 
^{wE{r)) (see FigI2J. Notice that hE ~ 6*"'^+'^'' grows at 
infinity - the QNM are not physical solutions themselves, 
but only approximate such solutions for a fixed r and for 



t ^ oo. The first 10 eigenvalues i?„ — — ijn are listed 
in Table I. 

Table I clearly indicates that ujn ^ 1 and 7„ for 
growing n. It seems that the QNM can be obtained for 
any n, thus comprising an infinite family. The values 
of LUn coincide well with the positions of the resonance 
peaks shown in Fig^ 

QNM — qualitative analysis.— The existence of the 
QNM of the BPS monopole can be qualitatively un- 
derstood as follows. Let us consider q introduced in 
Eqs.lO),® as a free parameter, q g [0, 1], and denote the 
corresponding solutions by w^g-j and For q — the 
equations decouple. Setting hny){t,r) — e^*"*(C+/i-(- -I- 
C-h ) with constant C±, Eq.l© is then solved by 

h± (r) = (coth r T iuj)e^''^''. (9) 

Eq.Q with iy(o)(<,r') = e^*'^*w(r) reduces to the eigen- 
value problem 

[-d^ + 72 jw{r)^^Mr)- (10) 

The potential in this equation has an attractive Coulom- 
bian tail, since it behaves as 1 — 2/r -|- 0{e~^) for 
r — > 00, thus there are infinitely many bound states, 
w{r) — Wn{r) for lu'^ = u)'^, n = 1,2,.... Several low 
lying LOnS are given in Table II. The n-th eigenfunc- 
tion Wn^r) has n — 1 nodes in the interval r g [0,oo) 
and can be normalized by the condition w'^dr — 1. 
For large n the t(;„(r)'s extend to the asymptotic region 
where the potential is Coulombian. As a result, they 
can be well approximated by solutions of the hydrogen 
atom problem. This implies that for n — > 00 one has 
u)'^ — 1 — l/n^ + 0{n^'^) and also that w„ ^ n~'^/^. More 
precisely, for w = 1 Ea. H10|l admits a limiting solution 
with infinitely many nodes, Wao{r), which itself is not 
normalizable, but Wooif) corresponds to the pointwise 
limit of 'Wn{r); i.e. for a fixed f 

lim Wn{r) — n^^^'^Wryo{r) \/r < r . (11) 

n — >oo 

Let us consider a solution of Eqs.jSJ,© given by 
/i(o) = and u'(o) = 5R(^„e~*""*'u;„(r)), where An is 
a constant. "Switching on" a small value of the coupling 
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TABLE I; The eigenvalues of the first ten QNM 



TABLE IL Values of and 7„ defined by Eqs. 10,(0 



f 



fO 



f 0.8473 0.5077 X fO~" 

2 0.9332 - iO.f 384 X fO"^ 

3 0.9637 - i 0.52f8 X fO-2 

4 0.9774 - i 0.2488 X fO~^ 

5 0.9846 -iO.f 375 X f 0-2 



7 

8 
9 

fO 



0.9888 
0.99f5 
0.9933 
0.9946 
0.9956 



i 0.8396 X fO-^ 
i 0.5499 X fO"^ 
i 0.3794 X fO"^ 
i0.273f X fO"^ 
i 0.2030 X fO"^ 



between the channels, (7 <C 1, the w-bound state will 
start loosing its energy to the /i-channel, where this en- 
ergy will be radiated to infinity. To find approximatively 
the corresponding solutions ft^) and W(q) by successive 
iterations, we solve first Eq.JoI) for by replacing w 
by 1/7(0) in its right hand side. The solution is regular at 
the origin and reduces to an outgoing wave at infinity: 
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Since there is now an outgoing flux of the /i-radiation, 
the energy of the w-bound state will be slowly decreasing. 
In the adiabatic approximation this process is described 
by a decrease of the amplitude of the w-field by replacing 
An — > An{t). To determine A„(t) we use the law of 
energy conservation, S^TJf = 0, whose integral form is 



d 
dt 



lim r'T^ . 



(12) 



Expanding Tjf up to terms quadratic in w and /i, the 
expression on the left in Ea. H12|) is proportional to -^A^ 
and determines the decrease of the bound state energy. 
The expression on the right is proportional to A^ and 
gives the energy flux at inflnity. Thus A„ = — 7„A„, 
from which An{t) = c„e~^"* with a constant c„. The 
coefficient 7„ is given by 
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(13) 



Summarizing, upon switching on a small value of q, the 
stationary bound state of the w-field, Wn{r), becomes 
quasistationary and is approximately given by 



w 



(9) 



-iE„t 



(14) 



where i?„ = LOn — ijm and hn{r) can be expressed by 
quadratures in terms of w„ (r) ■ Note the appearance of a 
complex energy. En, in Ea. ((T^ . Evaluating the integral 
in Ea. H13|l (see Table II) shows that 7„ decreases as n 
grows. This feature can be understood qualitatively by 
viewing the QNM as quasibound states of the massive 
w-field slowly decaying via an energy transfer to the h 
channel. For large n the coupling between the w and 
h channels becomes weaker and weaker, since it is pro- 
portional to WHwn ^ n^-^/^, and so the decay becomes 
more and more unlikely. Using the asymptotic relation 
(fTT|l in Eg. lfT^ yields 7„ ~ for n ^ 00. 

The approximative formulas derived above under the 
assumption q <ti I give in fact, when straightforwardly 



0.798 
0.057 



0.926 
0.010 



0.961 
0.0035 



0.984 
0.0009 



0.995 
0.0001 



extrapolated to q = 1, a good approximation for w„ — ijn 
(see Tables I, II). This approximation is actually getting 
better with growing n. We can therefore use the asymp- 
totic relations derived above to obtain for n ^ 1 



1 - ul =n + 0{n ^) , 7„ = bn' 



(15) 



Using the results of Table I, 6 w 0.2. 

Collective effect of the QNM.— The numerical sim- 
ulations in Ref . 4] of the non-linear temporal dynamics of 
the BPS monopole hit by a strong spherically symmet- 
ric pulse have shown that a considerable fraction of the 
energy received from the pulse is not radiated away im- 
mediately. It gets 'trapped' by the monopole and forms 
a long living quasiperiodic excitation that radiates very 
slowly, and whose late time behavior is to a large extent 
independent of the structure of the initial perturbation 
pulse, i.e., it shows certain type of universality. We can 
now offer an explanation to these observations. 

It is intuitively clear that, according to the decomposi- 
tion of the eigenfunctions of the linearized problem, there 
is a 'radiative' sector containing the massless /i-modes 
with Lu^ > and massive w- modes with uj^ > 1, and 
also a 'non-radiative' sector consisting of the w-modes 
with u!^ < 1. One expects that a part of the energy of 
the pulse received by the monopole will be distributed 
among the radiative modes and will be radiated away. 
However, as a generic perturbation will have a nonzero 
overlap also with the non-radiative modes, the remain- 
ing energy will get trapped in the non-radiative sector. 
Since it cannot be radiated away directly, a long living 
excitation on the monopole forms that will decay only 
due to a slow energy leakage to the radiative channels. 
In the terminology of black hole physics, the perturbed 
monopole keeps some of its 'hair' for a long time. 

Although initially the dynamic of the system will be 
nonlinear, one expects the nonlinear effects to become 
negligible when a sufficient amount of the received energy 
is radiated away. The linear approximation will therefore 
give a good description of the late time behavior, provid- 
ing also an explanation of the observed universality. 

Let t — he the starting point of linear regime, when 
the perturbed monopole is described by it;(0, r) = dW{r) 
and /i(0, r) = 5H{r). The subsequent temporal evolu- 
tion is determined to a large extent by the QNM, since 
they hold their energies for a long time. Using (|14|l . we 
therefore approximate the general solution for f > by 

w{t, r) = c„e-'""*-^"*«;„(r)^ , (16) 

and similarly for h{t,r). Here c„ = Wn{r) 6W{r)dr. 
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FIG. 3: The profile of w{t, 1) in Eg. 1161 corresponding to the 
initial step-like perturbation, SW{r) — d{l — r). 



This sum of damped oscillations corresponds to a long 
living excitation with the following properties. 

For a localized 6W the overlap coefficients c„ will be 
maximal for small n. Therefore terms with small n will 
dominate at first the sum in Ea. (|16|l . They will soon 
be damped, however, since their damping rates, 7„, are 
the largest. Terms with higher n will then become more 
important at later times. This 'dying out' of modes has 
indeed been observed in Ref. 4] (Fig. 6), and the values 
of frequencies uJn measured there are in good agreement 
with those given in Table I. 

For any t there is a number, k{t), determined by the 
condition ''^ 1, such that all terms with n <C fc in llt)|) 
are already damped, while those with n ^ k are not yet 
important, since their c„'s are small compared to c^. The 
sum will therefore be dominated by terms with n ^ k{t). 
Considering for simplicity the sum of only two of these 
terms, with frequencies Uk, (^k+i, gives beats with the 
base frequency Lu{t) — ^{uik+i + ojk) whose amplitude is 
modulated with the frequency r2(t) = ^{ojk+i—ojk)- This 
explains qualitatively the behavior shown in Fig|3| Since 
7fc decreases with k, the value of k{t) will grow with t, and 
so will do the base frequency uj{t). Since -fk ^ for 
large k, it follows that k{t) ~ t^^^ for large t. Using ((T3|l . 
we conclude that for large t one has 1 — to^(t) '-^ t^^^^, 
which explains the feature observed in Ref.Q (Fig. 5). In 
a similar way we obtain for the modulation frequency 

Even though each term in the sum Hlti|l tends to zero 
exponentially fast, the sum as a whole decreases consid- 
erably slower. Since for large t only terms with large 



n are relevant in (|16|l . we can use Eas. Hllll . (|15|) to ob- 
tain « 1 — l/2n^ and w„{r) « n^'^/^Waoir) for 
r < f — rw^dr ~ n^. For a localized SW this implies 
that c„ « Mnr'^l'^ with M ~ WaoSWdr. As a result, 
for t oo and for r < f{t) ^ n{t)'^ ^ t^^^ the solution 
(Uni reduces to w{t,r) = J\fwooirWiG(t)), with 

G(o = e- f; i»p(A-ii). (17, 

This shows that the late time dynamics is indeed uni- 
versal, since changing the initial conditions affects only 
the normalization J\f. The final task is to determine the 
asymptotic behavior of G{t). Transforming the sum H17|l 
to a contour integral and using the saddle point method, 
we find that \G{t)\ ~ t'^^^ for t oo. This t'^^^ expo- 
nent explains the feature observed in Ref.Q (Fig. 5). 

The perturbed monopole thus ends up in a long living 
breathing state dominated by the confined, slowly radi- 
ating massive modes of the gauge field. In the region 
r < t^/^ this state is characterized by modulated pulsa- 
tions whose base frequency approaches the vector boson 
mass, while the amplitude decreases as t~^/^ . 

The total energy of the breather can be obtained 
by summing over the QNM, E = c^w^e"^'*'"* ~ 
71^^ exp(— 26n^^i) ^ t^^/^. This includes the en- 
ergy of the confined massive w-modes and also that of 
the massless /i-radiation emitted by these modes. E de- 
creases, since there is a flux of the /i-radiation at infinity, 
S = E t~^/'^. This value agrees with the result of 
Specifically, Fig.2 in |1| shows the flux S ~ T'^/^ at fu- 
ture null inflnity, with T ^ t — r. By continuity, the flux 
through a 2-sphere of a very large but finite radius r is 
still S ~ since t is finite then, for t^r one 

has T ^ t, which agrees with out result S ^ t^^^^. 

In conclusion, we have found resonant excitations 
of the monopole and studied some of their applica- 
tions. Although we have considered above only the BPS 
monopole, the existence of similar resonances can also be 
shown in the case of a non-zero Higgs self-coupling A. 
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